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Abstract 

We consider the effects of temperature upon the superfluid phase of ultracold, weakly inter- 
acting bosons in a one dimensional optical lattice. We use a finite temperature treatment of 
the Bose-Hubbard model based upon the Hartree-Fock-Bogoliubov formalism, considering both 
a translationally invariant lattice and one with additional harmonic confinement. In both cases 
we observe an upward shift in the critical temperature for Bose condensation. For the case with 
additional harmonic confinement, this is in contrast with results for the uniform gas. 



I. INTRODUCTION 



Ultracold atoms confined within optical lattice potentials are of great current interest, 
both theo.etical.yll Q fl Q B a„d experimentally Q. 01 particula. interest are the 

connections to solid state regimes where optical lattices can be manufactured with extraor- 
dinary control so as to simulate more complicated, less perfect and, often, less manipulable 
systems, initially of interest in a condensed matter sphere (eg. Heisenberg or Ising Hamilto- 
nians )jay,[l^. Ultracold atoms in optical lattices have also been proposed as candidates 
for quantum information processing[ll|, which has applications in quantum cryptography 
and quantum computing. Various methods of loading Bose- Einstein concensates into optical 



lattices have been proposed 
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m 



15|, and this is now routinely performed jl^. as is 
the manipulation and control of the atoms in such a lattice jiol. 



The focus of this paper is on the microscopic treatment of Bose-Einstein condensates in a 
one dimensional optical lattice. We extend previous treatments at zero temperature BH to 
finite temperature using the Hartree-Fock-Bogoliubov mean-field treatment as applied to a 
discrete Bose-Hubbard model yielding modified Gross-Pitaevskii and Bogoliubov-de Gennes 
equations. This set of equations is solved for both the case of a translationally invariant 
lattice (no external trapping potential) and an inhomogeneous lattice (optical lattice in an 
external harmonic trapping potential). The use of the mean field based treatment means 
we are only considering the superfiuid phase for the atoms in the optical lattice. The model 
is not valid in the Mott Insulator regime. We use the model to estimate the superfiuid 
to normal phase transition temperature in each case. Thus one is able to obtain a phase 
diagram for the superfiuid and normal gas states in the low effective interaction strength 
limit. 



II. FORMALISM 



We consider a one-dimensional optical lattice with / lattice sites. We begin from the 
Bose-Hubbard Hamiltonian for atoms in a one-dimensional optical lattice 
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where J represents the couphng strength between adjacent lattice sites, V is the interaction 
potential acting between atoms on the same site, and di is the Bose field operator for 
the i^^ lattice site, and hi = d\ai . ei is the energy on each lattice site i due to the 
trapping potential. The usual commutation relations apply for the Bose field operator Sj. 
Assuming a macroscopic occupation of the ground state, we express the Bose annihilation 
operator for each lattice site i in terms of a complex mean field part Zi and a fiuctuation 
operator part 5i _17j, di = (^Zi + 5jj exp (^^), where /i is the eigenvalue for the generalised 
Gross-Pitaevskii equation to be discussed below, and take the self-consistent mean-field 
approximation such that SjSiSi = 2 (sjSi^ 5i + (^SiSi"^ Sj, etc. The condensate density is 

SjSi) is the excited population density. The anomalous 



given by ric = \zi\ and ru 
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The resulting Hamiltonian can now be diagonalised using the Bogoliubov (canonical) 
transformation 
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where and are the quasiparticle amphtudes, Uq are the quasiparticle excitation fre- 
quencies, and (dg) is the quasiparticle creation (annihilation) operator. This yields j^l a 
set of coupled equations comprising of a modified Gross-Pitaevskii equation 



fi'zi = e'iZi - {zi+i + Zi_i) + Kff {ric^Zi + 2niZi + rrtiZi) 



(4) 



and the Bogoliubov-de Gennes equations 



hu'^ul + dgZi = [2Vcs {ric^ + hi) - fi' + u'i\ u1 - + ul_^ - V^sz, 
-huj'^vl - d^z, = pKff (n., + h,) - + e[] vt - [vl, + vU] - V, 
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is the usual Bose distribution, the primed quantities are measured in units of J (which 
depends on the depth of the optical lattice), and where we have assumed without loss of 
generality the condensate amplitudes Zi to be real, /x is the energy eigenvalue for the modified 
Gross-Pitaevskii equation (and approximates the chemical potential closely for values of 
temperature well below the transition temperature) and huOq is the energy eigenvalue of the 
Bogoliubov-de Gennes equations. The c'^'s are necessary in order to ensure the orthogonality 
of the condensate with the excited states 0, Z is a fugacity term resulting from the 
difference between the true chemical potential /i-r, and the chemical potential as estimated 
using the eigenvalue /x corresponding to the ground state of the modified Gross-Pitaevskii 
equation. Thus Z = exp (/3 (/ly — /i)). 

In the case of a homogeneous gas with a large number of particles in the ground state 
(ie. 1), the fugacity may be approximated bvjisl 

Z = 1 + l/rie (8) 

In the lattice, this is not always a good approximation, since the condition rzc S> 1 does not 
always hold, and we will take Z = 1. This will lead to increased values for the excited atom 
population, and will therefore result in under-estimated values of the transition temperature. 

Superfiuid flow of the condensate occurs when there is a phase gradient. A phase gradient 
of the condensate modifles the hopping term of the Hamiltonian by the introduction of Peierls 
phase factors. The resulting energy shift may be estimated using second order perturbation 
theory. Since this energy shift is due entirely to the kinetic energy associated with the 
superfluid flow, and hence due to the superfluid fraction, it follows that the superfluid 
fraction fg may be calculated in terms of this energy shift 3| 

where J = iJJ2i=i (^^l+i^i ^ (^(^i+i^ T = JJ2i=i (^^I+i^^i + ^I^j+i j • The superfluid 
fraction can then be expressed in terms of the condensate and quasiparticle amplitudes 

/. = fi'^ - (10) 
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A. Translationally Invariant Lattice 



In the case of a translationally invariant lattice, periodic boundary conditions apply, and 
the quasi-particle amplitudes are given by 



(13) 



Furthermore the condensate amplitudes are equal for each site j and, since there is no 
trapping potential, = e = 0. Solving for the quasiparticle amplitudes we obtain 

Vz^ + AJsin^ (f) +huq 
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Thus from equations ()14|) . (fT3j) and (fTBj) . for the translationally invariant lattice, the super- 
fluid fraction is given by 



since q = 2Tij /la for 1 < J < / — 1- The condensate fraction for a given lattice site is given 
by fc = nd^o where no = N/I is the number of atoms per site. 

To calculate the condensate and superfiuid fraction, we first determine the condensate 
amplitude. Then defining 
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and 
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exp (^/5'e,(„_,) j - 1 

where j3' = J/ksT, ks is Boltzman's constant and where the subscript [n) refers to the 
variable in question at the n*'^ iteration, one can solve for \zf using an iterative scheme 
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As an initial guess, the value = y is used, and the calculations (fTH|) . (fT^ . (PU)) and 

(PT|l repeated until convergence is attained (ie. ~ l'^l(n-i) ^ Error Tolerance), or the 

maximum number of iterations is exceeded (divergent solution)). 

One first calculates gj(n), ej(„) and A''b£;^(„) using equations (fTH|) . (fT^ and The 
quasiparticle amplitudes (given by equations (fT^ and (fT3j) ) may be calculated using the 
equations 
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The condensate and superfluid fractions are then readily determined. 



B. Inhomogeneous Lattice 



The condensate amplitudes Zi are found by solving equation (jlj), where the trapping 
potential is given by = VL{i — (/ + l)/2)^ for site z, with VL = ^mu'^a^, and a is the 
inter-lattice spacing. The quasi-particle amplitudes for site i, w^and vfaie found by solving 
equations 

This set of equations can again be solved iteratively. In performing the calculation, we 
actually set the c'^'s to zero when solving the Bogoliubov-de Gennes equations, but do so 
in the Hartree-Fock basis, thus ensuring orthogonality of the ground state and the excited 
states 3]- The Hartree-Fock basis is given by the normalised solutions to the eigenvalue 
problem given by equation (jH), but where the ground state (zero energy solution) is excluded. 

First, let us rewrite the Bogoliubov-de Gennes equations ^ in matrix form 
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where 
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t = 2Kff {uc, +ni)-^l + n{i - (/ + l)/2f - J 



and 



(26) 
(27) 



Here JT" is defined as the operator acting on Zj, ul, and v\ as follows: 

+ (28) 

Now, let {zf jconstitute the eigenstates of equation (jH) with eigenvalues /i'^. We order 
the eigenvalues into ascending order, and order the normalised eigenstates accordingly, call 
these {^f }. The state Q corresponds to the Goldstone mode, we must exclude this in order 
to obtain the Hartree-Fock-Bogoliubov basis. Let us define the matrix 

U = 



where 



it 



(29) 



(30) 



Let v^^^jbe the solution to the Bogoliubov-de Gennes equations in the Hartree- 



Fock-Bogoliubov basis, then the solution to the matrix equation 
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U^CU U^MU 
-U^MU -U^tU 

gives the Bogoliubov quasiparticle amplitudes in the Hartree-Fock-Bogoliubov basis. To 
obtain the Bogoliubov quasiparticle amplitudes{M^, } one transforms back using 

i9 



^ "-HFB 
U^HFB 
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The iterative solution of these equations is then obtained numerically. 



III. RESULTS 

We present results for a lattice of depth of approximately IG.SEn, where the recoil energy 
is defined as En = h'^/Sma'^ with a the inter-site spacing. This is defined for the optical 



lattice parameters defined in reference |2fl| and using the band structure calculations of 



2l| 



A. Translationally Invariant Lattice 

In figures n and 121 we present results for a translationally invariant lattice with one atom 
and ten atoms per site respectively for various values of the effective interaction potential 
Vcs- In both cases ten lattice sites, with periodic boundary conditions, were used. One 
observes, a decrease in both the condensate and superfiuid fractions with temperature and 
we interpret the point at which these densities approach zero as indicative of the critical 
temperature Tc for the superfiuid to normal gas phase transition. This interpretation is 
supported by an examination of the low lying excitation spectrum which, for the case Ves = 
20, is shown in panel (b) of figure |21 The "softening" of the modes, indicative of the phase 
transition, is clearly seen and coincides with the transition temperature obtained by noting 
the temperature at which the condensate and superfiuid densities approach zero. 

We obtain such a transition temperature for each value of the effective interaction poten- 
tial Ves, producing a phase diagram as shown in figures [H^c) and|2l^c), where the superfiuid 
phase and the normal phase are as indicated, the superfiuid lying to the left of the curve. 
The transition to the Mott insulator phase cannot be determined by this analysis, but one 
would expect the transition from the superfiuid phase to the Mott insulator phase to occur 
when the site coupling strength J is decreased below a certain point, depending on the 
on-site interaction strength Ves- In practice the Mott insulator phase transition will occur 
when Ves = V/J exceeds some critical value. Thus a transition to the Mott insulator phase 
is also possible by increasing the on-site interaction strength V for a given coupling strength 
J. We note here that, in all instances, the transition temperature increases with the effec- 
tive interaction potential, as is indeed the case with a homogeneous Bose gas (ie. a Bose 



gas in the absence of an optical lattice and of a confining potential) 



This reentrant 
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FIG. 1: Translationally invariant lattice with ten atoms and ten sites. Overall condensate (a) and 
superfluid (b) fractions as a function of temperature. The corresponding phase diagram is shown 
in panel (c). 

behaviour has also been predicted for the translationally invariant lattice by Kleinert et al. 



2J| . Our results are consistent with their conclusions. In addition, they predict a reduction 
of the critical temperature as the effective interaction strength is increased further. We, 
however, are unable to explore this regime as it extends beyond the validity of our model. 
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FIG. 2: Translationally invariant lattice with ten atoms per site. Condensate Fraction (a), and 
excitation spectrum as a function of temperature (b) with the corresponding phase diagram (c). 

B. Inhomogeneous Lattice 



In this section we present results for the case of an optical lattice in the presence of a 
harmonic trapping potential. The condensate fraction fc and the superfluid fraction fg are 
evaluated as a function of the effective potential Ves and the temperature. 

Figure El shows plots of the number of condensate atoms, excited atoms and of the total 
number of atoms for each lattice site for the case of an inhomogeneous optical lattice con- 
sisting of forty one lattice sites (odd case) with ten atoms in total, for Ves = 1 at various 
temperatures ranging from T = nK to T = 1.6 nK. At zero temperature, the condensate 
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FIG. 3: Inhomogeneous lattice with forty one lattice sites and ten atoms in total for V^fi = 1, at 

(a) T = nK, (b) T = 1 nK (c) T = 1.6 nK. Circles represent the number of condensate atoms, 
squares the number of excited atoms and crosses the total number of atoms 

atom distribution is bell-shaped, peaked at the central lattice site. There is a small quantum 
depletion even at zero temperature, and the distribution of excited atoms is shaped as a 
bimodal distribution, centred about the central lattice site. As the temperature increases, 
the condensate population decreases, but the distribution still remains bell-shaped, and the 
excited population increases. In figure IH we present the overall condensate and superfiuid 
fractions and the corresponding phase diagram for optical lattices in a harmonic potential 
consisting of forty lattice sites (even case). Panel (a) corresponds to the condensate fraction, 

(b) the superfiuid fraction, and (c) the phase diagram. We note that for higher tempera- 
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FIG. 4: Overall condensate and superfluid fractions as a function of temperature, and the corre- 
sponding phase diagram for an optical lattice in a harmonic potential with forty lattice sites (even 
case) and ten atoms, (a) corresponds to the condensate fraction , (b) to the superfluid fraction 
and (c) to the phase diagram. 



tures, calculations performed using forty one lattice sites begin to show a marked difference. 
This is indicative of the fact that we are pushing the bounds of validity of our model. In 
particular we are seeing significant finite size effects, which affect the value of the chemical 
potential and, ultimately, a failure of the mean field approximation. The dotted continuation 
of the lines in figure 0] indicate where our calculations become unreliable, but are included 
for completness. We are still able to use our model to obtain an estimate of the critical tem- 
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perature and hence the trend in its dependence upon the effective interaction strength. We 
therefore conclude from figured that the superfiuid to normal phase transition temperature 
increases with increasing Ves- It is clear, then, that the shift in critical temperature with 
effective interaction potential is positive definite for a Bose gas in a one-dimensional optical 
lattice, regardless of whether the system is confined in a (harmonic) trapping potential or 
not. This is in contrast to the case of a three-dimensional Bose gas, where ATc changes sign 
for the trapped gas. 



IV. CONCLUSIONS 



We have applied the descretized Hartree-Fock-Bogoliubov formulation to the Bose- 
Hubbard model in order to calculate the dependence of the condensate and superfiuid frac- 
tions on the temperature. We have used this to estimate the critical temperature for the 
superfiuid to normal phase transition for both a translationally invariant optical lattice (no 
external trap present), and an inhomogeneous optical lattice (contained within an external 
harmonic trap). This has enabled us to investigate the phase diagram for both cases and 
we observe that the transition temperature increases with increasing effective interaction 
potential V^g. Unlike the homogeneous case with no optical lattice, this positive shift in 
the critical temperature with interaction strength is present in both the translationally in- 
variant case and when a parabolic confining potential is imposed. In the homogeneos gas 
the shift in the critical temperature is only positive in the absence of a confining poten- 
tial. These conclusions are consistent with previous work for the translationally invariant 
case ■ extending this result to include parabolic confinement. 
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